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Abstract 



We consider the Affleck-Dine mechanism for leptogenesis in the minimal MSSM 
with Dirac or Pseudo-Dirac neutrinos. The rolling of scalars along D-flat di- 
rections generates a left-right asymmetry in the sneutrino sector, only the left 
part of which is transferred to a baryon asymmetry via sphaleron transitions. 
C*~> ' In the pure Dirac case the baryon asymmetry of the Universe is thus mirrored 

by an equal and opposite asymmetry in the leptons. The mechanism is also 
^-j- \ found to work when the neutrinos are pseudo-Dirac. No additional field needs 

I ' to be added to the MSSM other than the right-handed neutrino. 

o 

1 Introduction 

q h i It was noticed some time ago that leptogenesis, far from requiring lepton- 

^ ■ number violating Majorana masses as in the original scenario [1], can in fact 

<D \ be implemented in the SM with purely Dirac neutrinos [2]. This mechanism, 

' called neutrinogenesis, relies on the fact that (i?+L)-violating transitions leave 

■ the right-handed sector unaffected [3]; as long as left-right equilibrating pro- 

^ , cesses are small enough to be out of equilibrium, it is possible to 'hide' a right- 

!-h ' handed lepton asymmetry from the sphaleron transitions. The idea of hiding 

lepton number in inert species has a long history [4] but works particularly 
effectively for the neutrinos; indeed it was shown in [2] that Dirac neutrinos 
easily satisfy this condition. The temporary left-handed lepton asymmetry 
can thus be processed before the electroweak phase transition into today's 
observed baryon asymmetry. It is only well after the phase transition that 
the neutrinos' Yukawa couplings come into equilibrium, by which time the 
sphalerons are quenched and the baryon asymmetry is locked in. The useful- 
ness of this idea lies in the fact that Dirac neutrinos of the right size can arise 
in models where GUT scale degrees of freedom are integrated out because 
Yukawa couplings in such models can be naturally suppressed by factors of 
M\y /Mqut- This possibility has received increased interest recently in the 
context of supergravity and effective models of string theories [5-10]. The fact 
that baryogenesis is also possible then leaves open the intriguing possibility 
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that B — L is conserved in Nature or that neutrinos are in fact pseudo-Dirac 
rather than Majorana. 

However the toy model used in [2] used an additional heavy Higgs-like 
doublet, because the scenario worked by 'drift and decay' as in original lep- 
togenesis. In the present paper we point out that the Affleck-Dine (AD) 
mechanism [11] allows an extremely efficient implemententation of neutrino- 
genesis in just the MSSM with Dirac neutrinos. A {pL — vr) current can be 
produced through the rolling of scalars along their D-Rat directions; although 
lepton-number is conserved, only left-handed lepton number can be converted 
to a baryon number through sphalerons, and the right-handed component is 
hidden by the smallness of the Yukawa coupling as before. The B + L number 
of the Universe is thus mirrored by an equal and opposite right-handed lepton 
number, until the right-handed (s)neutrino oscillations decay long after the 
electroweak phase transition. We should mention that AD neutrinogenesis 
was proposed in ref. [12]. However in that work the AD field was considered 
to be an additional scalar field that was either Higgs-like, with SU(2) number, 
or a singlet appearing in higher order non-renormalizable interactions. The 
implementation here using only the D-flat directions of the MSSM itself can 
be thought of as the minimal realisation of AD neutrinogenesis in the context 
of supersymmetry. Apart from testing whether AD can work in a B — L pre- 
serving MSSM this minimal scenario is naturally very predictive, because for 
example the CP violating potential arises from soft-supersymmetry breaking 
trilinear terms, the coupling of the AD field to the neutrinos is given by the 
neutrino mass, and so on. 

We first discuss the evolution equations for the left-right (L — R) asym- 
metry as a result of the Dirac mass term added to the MSSM superpotential. 
We will then evolve numerically the L — R asymmetry; using the equilibrium 
relations between B and L, this asymmetry will be converted to a baryon 
asymmetry. Finally we discuss the pseudo-Dirac and 'weak see-saw' cases. If 
neutrinogenesis is to work a bound on the Yukawa couplings results from the 
requirement that the oscillations of the AD field decay after the electroweak 
phase transition. The nett result is a constraint on the size of the additional 
Majorana mass which turns out to be 



This bound is the strongest, being slightly more severe than the requirement 
that the baryon number generated be large enough. In short this means 
that the AD neutrinogenesis is able to operate for all 'reasonable' Dirac and 
pseudo-Dirac scenarios in the neutrino mass sector. It even works for mild 
see-saw cases although the latter are probably excluded by nucleosynthesis. 



Let us first introduce the right-handed neutrino superfield, N, and add a Dirac 
mass term for the neutrinos in the superpotential: 




2 The Superpotential 



V s D AI/e^HVN , 



(1) 
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where L is the left-handed lepton doublet and H u is the up-type Higgs. The 
gauge invariants LH U and N will be the important D-flat directions [13]; let 
us parameterise them as 
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where non-bold letters stand for the scalar part of the superfields. The (SUSY- 
conserving) scalar potential arising from the added Dirac mass term is 
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|<A 2 | 2 + |A| 2 |P<5 



dW 



dH b 
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dW 



(3) 
(4) 



These tiny F-term contributions lift the LH U and N flat directions very 
slightly In usual AD the directions considered would be both D and F- 
flat at the renormalizable level, and the flat direction would be lifted only 
by non-renormalizable (i.e. higher dimension) operators. Here the directions 
would be F-flat as well but for the neutrino mass. It is only because of the 
smallness of the latter that we can hope to send the field out to large enough 
VEVs to generate asymmetries. Note that conversely when we go on later to 
consider pseudo-Dirac neutrinos, then the scenario begins to run into difficulty 
if there is any significant see-saw effect at work in the neutrino mass matrices: 
a significant see-saw would imply larger Dirac Yukawa couplings and lift these 
directions more. 

The AD mechanism of course requires additional CP violation, and here 
it comes from the soft-breaking sector; 



V SB 



\ 2 + ml 



v\ + (\a(p 2 v + /i.e.) 



(5) 



As was discussed in [13,14], soft-breaking terms also get a contribution from 
the non-zero Hubble constant in the early Universe and this is crucial as it 
drives the fields out to large values during inflation. We parameterise these 
as 

V H = -c^H 2 - c v \vf + (Xc H H(P 2 9 + h.c.) . (6) 
The overall potential for the scalar fields is thus 



V = V f + Vsb + Vh 

= {ml - H H 2 ) + (ml - c v H 2 ) \l\ 2 + (A(a + c H H)<j) 2 D + h.c.) 

+ ^V| 2 + |A| 2 |P0| 2 . (7) 



For the flat directions to develop large expectations values during inflation, 
we need at least one of the fields to have a negative effective mass squared; 
here we will consider 

(m 2 - C4 ,H 2 ) < . (8) 
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Thus the Hubble induced terms in eq.© push the fields far from the origin. 
They are also important in introducing a time-dependence into the potential 
which guarantees that the AD mechanism will be in operation. To see this 
note that without these terms the potential always has a minimum where 
sign(Aac/) 2 v) = — 1 as this is the only trilinear term. Without an initial kick 
the fields will simply roll down this valley and no nett lepton currents of any 
kind will be generated. The Hubble induced terms mean that the minimum is 
now at sign(A(a + chH)4> 2 v) = —1. Thus even if ch is real the phases of the 
fields have to become time dependent to track the instantaneous minimum: in 
effect the Hubble constant should kick the AD field for us wherever it starts 
out. 



3 The dynamics of the asymmetry 

To see these effects we proceed to examine how the asymmetry develops. First 
write the lepton number til as a sum of its right-handed and left-handed parts: 

n L = nf } + nf ] (9) 

with n L and n L being in terms of our scalar fields 

(L) i 
L 2 

= v — u*i^j . (10) 

The evolution equation for is: 

* + 3ff*+|£ = (11) 
and analogously for v. Now using ea. (|l(J|) in ea. (|lH) and its conjugate, we find 



4 L) + 3Hn[ L) = Im I ^ I , (12) 

and again analogously for v. From eq.0 we see that the only imaginary 
terms are the a-terms and hence 



dV 



n 1 ^' + 3Hn\ ' = 21m (\a(f) 2 i>) 
njf 5 + 3Hn^ } = -21m (Xacf) 2 ^) . (13) 

We can see that lepton number eq.@ is conserved, 

h L + mn L = j t [n^ ] + nf } ) + 3H {n^ ] + nf ] ) = 0, (14) 

but that the left-right asymmetry, — = ulr, has a non-trivial evolu- 
tion: 

h LR + 3Hn LR = Alm(\a(j) 2 D) . (15) 
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The hope then is that this time dependence and CP asymmetry in the 
potential will generate a nett ulr. If it does so this will feed through to the 
baryons via sphalerons. Before continuing we briefly establish the relation 
between n L and the baryon number (before taking account of the effect 

( Ft} 

of sphalerons, n L is initially half the value of ulr generated during the 
evolution of ea. (|15|l ): the equilibrium ratio between lepton and baryon number 
under rapid sphaleron transitions was calculated in ref. [15,16] for an SM like 
structure. In the present case we have an out-of-equilibrium right-handed 
Dirac neutrino in the analysis which simply holds a nett B — L and remains 

( R} 

completely inert; therefore we can set (B — L) = n L where L is the sum of 
all the leptons in equilibrium (i.e. excluding the right handed neutrino). In 
the MSSM we have an additional charged Higgs which changes the result from 
the SM; repeating the chemical potential analysis and assigning a chemical 
potential ^ B _i we find that above the electroweak phase transition for m 
Higgs fields 

Y = IQ^b-l + ( 20 + 2m)/iy 
B-L = 26/i B _£ + 16/iy (16) 

while below it 

Q = 16/i B _£ + (44 + Am)/J, Y 
B-L = 26/^ + 16^. (17) 

Imposing (B — L) = n\ this then translates into the following ratios: 

R j _ 4(6+m) (R) _ _8_ (R) rp rp 

D — 1j — 66+13m >l L — 23 11 L 1 ^ 1 ' 

4(9+m) (R) 44_ (1 

Ul+13m n L ~ 137 " L 



66+13m"L — 23' L x ^ f-,o\ 

R r 4(9+m) (R) _Uj R ) T<-T 1 ' 



where L is the total lepton number including the right handed neutrinos and 
m = 2 in the MSSM. (The only effect of the charged Higgs of the MSSM is 
to change the denominator of the last expression to 111 + 13m rather than 
98 + 13m.) Note that in this pure Dirac case B — L is conserved and the final 
baryon number is approximately the LSP density, given by udm = f «B) 
reminiscent of ideas pursued in a number of works [17-27]. This suggests the 
right handed sneutrino as the preferred LSP candidate since the LSP mass 
would then have to be of order lGeV: 

8 &.DM rin x 
m DM = -^—Q^~ mb ■ ( 9 ) 

Returning now to the dynamical evolution first note that the forcing term 
for Ulr evolution is time dependent and this is true generically when H ^> 
m 3/2 even if the fields are sitting in the minimum because the minimum of 
the potential depends on H(t). Indeed we have to minimise 

V = (ml -c^l^ + iml-c.H 2 )^ 2 -2\X(a + c H H)\\^\ 2 \D\ 
+ |A| 2 U 2 | 2 + |A| 2 |^| 2 . (20) 
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Taking the coefficient of \u\ positive, and for |c„| H 2 3> m 2 and \c^\ H 2 S> rrA 
the minimum of the potential is given by 



One obvious difference between the present case and the more usual one in [13] 
is that here the initial field values have a lower bound of order </> m i n — a/A ~ 
10 GeV even if H itself is much smaller than this value. 

Let us now sketch the evolution from early after inflation to the electro- 
weak phase transition; we will confirm the picture with a numerical solution 
of the equations of motion as we go along. The numerically evolving ulr is 
shown in fig.(l). 

The mechanism requires that initially the fields are drawn far along the 
flat directions during inflation. This in turn requires the coefficient of \4>\ 2 
in eq.Q to be negative as we have seen above. Inflation is then followed by 
an era of inflaton oscillation, during which the Universe is matter-dominated 
{H ~ 2/3t). At this early stage the time dependence of the Hubble constant is 
still important in the evolution of the fields. Indeed both (f> and v are following 
their evolution equations and, since H 3> m^,m=, their motion is dominated 
by the falling value of the H 2 Hubble induced mass-squared terms. 

The behaviour of the fields in this phase is a major difference between 
the scenario we are considering here and the original AD scenario. The flat 
direction here is lifted by renormalizable terms; it is easy to see from the 
equations of motion that the distance of the fields from the instantaneous 
minimum drops as t ; but eq. (|21j) tells us that the minimum itself drops as 
t^ 1 as well, so that during this matter dominated phase the fields are relatively 
undamped. We expect to see a long transient period and as our numerical 
analysis shows, this is indeed the case. To get a crude understanding of 
the behaviour in this phase, we can estimate the maximum amplitudes of 
the fields by assuming that the energy is constant in a co-moving volume: 
i? 3 F 2 (/4 ax = const. This gives </> max = const which in turn gives = const. 
The detailed behaviour of niR is still rather complicated at this point; the 
fields are not yet executing regular cycles, and the current ulr is rapidly 
flipping sign as the fields change their sense of rotation around the origin. 
This behaviour is in agreement with the arguments of refs. [13] where it was 
noted that only with nonrenormalizable terms of dimension 4 or higher do the 
AD fields follow the instantaneous mininum closely. 

It is during matter-domination that the important H ~ w.3/2 mark is 
reached. Below this point the Hubble induced terms in the effective potential 
become irrelevant to the evolution which is now dominated by the mass terms, 
and the behaviour changes markedly. The time dependence of the fields be- 
comes more amenable to analytic approximaton now since the equations of 
motion are nearly linear; indeed we may immediately use constancy of energy 
in a comoving volume argument to infer that .R 3 m 2 </> max = const where m is 




(21) 




chH <C |a| , 
chH » |a| . 
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the mass of whichever field we are considering. If H = b/t this then suggests 



36 

--, (22) 



which we indeed confirm numerically. In matter domination b = 2/3 so that 
riLR drops as t -1 . This can clearly be seen in fig.(l). Slightly later, but before 
reheating, the fields begin to exhibit "canonical" AD behaviour where they 
do execute regular cycles. Here we may approximate the real and imaginary 
components of fields as t k sin(mt). Linearizing the equations of motion we 
find 

k(k - l)sin(m^t) Ikm^ cos(m^t) Sbm^ cos(m^) _ 

P + t + t " U ' [ 6) 

and then neglecting 1/i 2 terms we find k = —36/2 agreeing with the above. 

The fields stay in this phase until reheating when the Universe becomes 
radiation dominated. We assume that reheating happens for Tr = 10 9 GeV 
when the Hubble constant is H ~ T|/M P/ ~ 1 GeV. We now have H = l/2i, 
which gives k = —3/4, in accord with the late time behaviour of the Bessel 
function solutions of ref. [11]. The current then drops as ulr ~ i -3 / 2 as 
regular matter until the electroweak phase transition. This behaviour is again 
confirmed by the numerical treatment: in particular fig.(l) shows the initial 
transients, the switch to AD behaviour, the t _1 and the i~ 3 / 2 behaviour. 



It is in this final phase as the fields are rolling down to the minimum, 
that they capture a left-right asymmetry that is constant relative to the en- 
tropy. A positive ulr means for instance that there is instantaneously more 
left-handed sneutrinos and right-handed anti-sneutrinos than left-handed anti- 
sneutrinos and right-handed sneutrinos, respectively: the left-handed sneutri- 
nos are quickly turned into left-handed neutrinos through gaugino interac- 
tions. This can either go by decay with V ~ g^m- or at high temperatures by 
scattering whose rate is 

T f 2 - — T 5 (24) 

m w,B 

where the masses are understood to be thermal ones. All of the contributions 
are of the same order during the period we are considering when T ~ My/ 
and so sneutrino<-^neutrino conversion is in equilibrium. The sphaleron tran- 
sitions transfer the left-handed neutrino asymmetry into a baryon asymmetry 
as described above. Above the electroweak phase transition this happens on 
a timescale of order TeV -1 which is essentially instantaneous; after the elec- 
troweak transition the sphalerons are switched off and the non-zero baryon 
number is frozen in [3,28,29]. Throughout, the right-handed (s)neutrinos 
remain inert until their oscillations decay through the coupling to the neu- 
tralino. We should therefore ensure that this happens at a temperature much 
lower than the electroweak transition temperature T ew : the life-time of the 
sneutrino decay is 

47T 1 

^ " A 2 rn~^ Higgs (25) 
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log (n LR ) 




Figure 1: Time evolution of the generated LR asymmetry. Parameters and initial 
conditions are as follows: = 600 GeV, m~ = 500 GeV, a = e a6i 100 GeV, = 1, 

Cl = 0.8, C H = 0, A = 10~ 12 , (f)(t in ) = l |0| min (t^ "(tin) = \u\ min (t in ), 0=^ = 0, 

where the minima are given by the expressions in the text. The added line is matter 
evolution during radiation domination, £~ 3 / 2 . The behaviour of the (J) field is also 
shown for early (shortly before H ~ lOOGeV) and late (post-reheating) times. 
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where -Bffiggs < 1 is the fraction of the neutralino that is made of the up 
Higgs. To be conservative we take B = 1, and find that t$ r > H^ 1 for all 
T > (j^n) 100 MeV. Note for later use that the mechanism stops working 
when 

A > 10- 9 (26) 

because the oscillations are damped before the electroweak phase transition 
takes place. 

The discsussion above of course assumes that the LSP is the usual neu- 
tralino. In the possibility we've mentionned above that the right-handed sneu- 
trino itself be the LSP, the decay time for the sneutrino oscillations is even 
later. 



4 The baryon asymmetry 

Having established the dynamical behaviour of the fields in some detail, let 
us now turn to the requirements to successfully generate a baryon number. 
First we have seen that the oscillations in tilr remain constant for an initial 
transient period when the inflaton oscillations are scaling like matter. They 
begin to scale like matter as well once the Hubble constant reaches ~ ^3/2- 
Therefore in order to evaluate the eventual tilr asymmetry it is most useful 
to consider the relative densities when H ~ m-3/2) the density in coherent 
inflaton oscillations is pi ~ H 2 Mp ~ my 2 Mp. At this stage the field VEVs 
are of order <fi, v ~ |a/A| as in ea. (|21|) . so that the energy density in their 
oscillations is of order p^^ ~ mg^ 2 |a/A| 2 . Since the latter also behaves like 
regular matter, we can use it to keep track of nm until the time of reheating: 



\a/\\ 



Pi M 



2 ■ (27) 
p 



From reheating onwards it is the ratio with entropy that remains constant. 
Since = m^n^ - that is given by 



|q/A| 2 T R 
Ml m 



lO" 9 



lOOGeV 



10 



-12 



5 The pseudo-Dirac and mild see-saw cases 

It is interesting to implement this scenario in the more general case where 
Majorana mass terms are included in the superpotential 

WDNAL a yH u t + M Ji NN + -^(LH 11 ) 2 . (29) 

Such additional terms can arise in the same manner as the Dirac terms in 
the supergravity scenarios considered in ref. [5-10]; essentially the pure Dirac 
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models require symmetries to prevent Major ana masses for the right-handed 
neutrinos that can be relaxed to allow non-renormalizable operators such as 
HuHdNN/Mct/T- For example such an operator could lead to left-handed 
Majorana masses Ml ~ 3 x 10~ 5 — 7 x 10~ 4 eV in the models considered 
in ref. [10]. In order to present as general a discussion as possible we will 
consider Ml,r to be arbitrary parameters and consider the question of when 
AD neutrinogenesis can work. We will mainly focus on Mr since Ml > Xv 
would give 6 active neutrinos which is certainly ruled out by nucleosynthesis; 
other than this we will consider Ml,r to be free parameters. Throughout the 
following discussion we shall assume that mass-squared differences given by 
measured neutrino oscillations are indicative of the actual masses. 

The most immediate concern is how the new terms could affect the classical 
dynamics. Assuming for the moment that Ml = 0, the superpotential leads 
to the scalar potential (J7J): 

V = (m$-c*# 2 )M 2 + (ml- CN H 2 + AM 2 R )\~v\ 2 
+ (X(a + c H H)(p 2 v + XM R( j) 2 u* + h.c.) 

+ ^ 2 | 2 + |A|W . (30) 

Clearly a new trilinear term has appeared which, following ea. ()ll|) . could affect 
the dynamics of the fields, and thus of the asymmetry if Mr > a. However 
as we shall see such large values are not relevant for the AD scenario here. 

Next, the lepton- number violating interactions introduced by the Majo- 
rana mass have to be constrained such that they do not erase the asymmetry. 
The rate of these interactions when Mr ^ is given by an exchange with vr 
which for T 3> Mr is 



A 4 M 2 



LV — 



(31) 



T ' 

Demanding that this rate be smaller than H for the duration of neutrinogen- 
esis imposes 

A a Mj|<|k (32) 

which taking T ew = lOOGeV gives 

X 2 M R < 3 x 10~ 7 GeV. (33) 

As we shall see the resulting bounds are not very constraining. For Ml the 
lepton number violating exchanges are now suppressed only by the gauge 
couplings, c/|, rather than A 4 ; however this still gives an uninteresting bound 
M L < 0.3MeV. 

A different and more constraining bound comes from the Dirac Yukawas 
themselves; again considering Ml = 0, the light neutrino mass, given by 



'Ml + (2Xv) 2 - Mr 
m v = ^ , (34) 

can require a A substantially larger than 10 -12 due to see-saw effects once 
the Majorana contributions become dominant; any bound on A coupled with 
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the estimated neutrino mass puts an indirect bound on Mr. One such bound 
comes from the fact that if A > 10~ 8 then left- and right-handed neutrino can 
equilibrate above the electroweak phase transition, destroying any left-right 
asymmetry [2]. However as we saw a stricter bound, 



A < lO" 9 



(35) 



comes from the requirement that the oscillations remain undamped until after 
the electroweak phase transition. In addition we of course require that the 
produced baryon number is sufficient; from ea. ()28j) we see that this is least 
constraining for maximum reheat temperature. Assuming that the gravitino 
reheat bound is satisfied Tr < 10 9 GeV and that ur/s ~ 10~ 10 we find the 
bound is (coincidentally) almost the same as in eq. (|33|). A < 3 x 1(T 9 . The 
added Majorana mass Mr must therefore respect both the constraints in (|33j) 
and (|35|) . It is the latter which is most constraining; it gives 



This bound then prevents the see-saw mechanism from operating in all its 
glory, but a weak see-saw effect may still be present in the neutrinos while still 
allowing AD neutrinogenesis to work with these flat directions. In addition 
note that the bound means that the evolution is virtually unnaffected by the 
presence of of the Majorana terms since Mr is indeed much smaller than both 
H and the trilinear term a driving the dynamics. 

In this context, following eq. (|28j) . there are two options for explaining why 
ur/s ~ 10 -10 . The first is the possibility that Mr ~ lMeV and that there is 
a weak see-saw mechanism in operation. This seems rather unnatural since it 
introduces an additional mass-scale that itself requires explanation. Moreover, 
it has been argued that a lMeV sterile neutrino with such a (relatively) large 
mixing ang le (sin 2 29 = 2 x 1(T 7 ) would cause large amounts of sterile neutrino 
dark matter to be produced [30]. This scenario would then be forbidden by 
overclosure (and possibly other cosmological constraints - see [30]). In such 
a case, the remaining possibility is that the reheat temperature was of order 
ITeV 1 . 

6 Conclusion 

In this paper we presented a minimal version of neutrinogenesis with Dirac 
sneutrinos in the MSSM, and showed that it can generate the observed baryon 
asymmetry of the Universe. No new fields need to be added to the MSSM 
apart from right-handed neutrinos. The mechanism works by first generating a 
V L — v r asymmetry using the AD mechanism, with Z)-flat directions involving 
sneutrinos and Higgses playing the roll of the Affleck-Dine fields. The flat 

: It has been noted that such a low reheating temperature could arise in supersymmetry [31], if 
one of the D and F-flat directions acquires a VEV of order the Planck scale. (The flat directions 
required for the mechanism here do not get large enough VEVs to affect the reheat temperature 
by more than an order of magnitude) . 




(36) 
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directions are appropriately lifted during inflation by the inclusion of finite- 
energy density SUSY-breaking terms which drives the VEVs to large values. 
As long as left-right equilibration is out of equilibrium before the electroweak 
phase transition (resulting in a bound on the Dirac neutrino Yukawas), the 
nett left-handed lepton number can drive sphaleron transitions and ultimately 
create the observed baryon asymmetry. We also showed that the conditions 
on the smallness of the Yukawa couplings still allows the mechanism to be 
implemented for pseudo-Dirac neutrinos, and can in fact support a weak see- 
saw mechanism. 
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